Math 155: Homework # 9, due on Friday, March 13 (6 questions) 

[1] Consider a linear, position-invariant image degradation system with impulse response 

h(x -a,y-P) = e -[(*-°) a +(»-0) 2 ]. 


Suppose that the input to the system is an image consisting of a line of infinitesimal width located 
at x = a, and modeled by f(x,y) = 5(x — a), where 5 is the impulse function. Assuming no noise, 
what is the output image g(x, y) ? 

Hint: Use equation (5.5-13) and also that 
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where y and er are constants. 

[2] (a) Consider the radial function ho(r) = e _r "/ 2 °' 2 , where r 2 = x 2 + y 2 . Compute the 2D 
Laplacian V 2 [/io(\/x 2 + y 2 )]. 

(b) A certain X-ray imaging geometry produces a blurring degradation that can be modeled as 
the convolution of the sensed image with spatial, circularly symmetric function 


h(r)= (r 2 — 2a 2 )/a 4 


-r 2 /2o 


where r 2 = x 2 + y 2 . Show that the degradation in the frequency domain is given by the expression 

H(u,v) = -8vr ^a 2 {u 2 +v 2 )e~ 2n2(j2 ^ 2W) . 


[3] Parametric Wiener Filter 

(a) Implement a motion blurring filter as in problem [6], Hw # 8. 

(b) Blur image 5.26(a) in the +45o direction using T = 1, as in Fig. 5.26(b) (a = b = 0.1). 

(c) Add a small amount of Gaussian noise of 0 mean to the blurred image. For this, you can 
use the Matlab command: 

J = imnoise(l,’gaussian’,0,v) adds Gaussian white noise of mean 0 and variance v to the image 
I (the parameter v needs to be specified). The default is zero mean noise with 0.01 variance. 

(d) Restore the image using the parametric Wiener filter given by 


F(u, v) 


- 1 \H(u,v)\ 2 - 

-H(u,v) \H(u,v)\ 2 + K- 


G(u,v), 


where K is a specified constant, chosen to obtain best visual results. 

(note that you may have to avoid division by zero in the above formula). 


[4] (a) Recall the definition of the convolution f*g(x, y) in continuous variables and two dimensions, 
(b) Show that V 2 (/ * g) = f * (V 2 g), where V 2 denotes the Laplace operator in (x,y). 

[5] (a) Read Section 10.2 (pages 692-717). 

(b) Give the main steps of edge finding using the zero-crossings of the Laplacian of Gaussian 
(read pages 714-717). 

(c) Implement and apply this method to the angiogram image from Fig. 10.15(a) (as in Fig. 
10.22 for the house image, but you will apply your method to the angiogram image). 
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[6] Recall the ID Laplacian of Gaussian (LoG) operator 


V 2 /r(r) 




(a) Show that the average value of the LoG operator V 2 /r is zero. Hint: use the identities 
- 0=5 /+£ e~ r2 /^ 2 Ur = 1, a 2 = -X=* /+°° r 2 e~ r2 /^dr. 
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(b) Prove that the average value of any function convolved with this operator is also zero. Hint: 
use the frequency domain. 
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